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We discuss the sign problem in QCD at nonzero chemical potential and 
its relation with chiral symmetry breaking and the spectrum of the Dirac 
operator using the framework of chiral random matrix theory. We show 
that the Banks-Casher formula is not valid for theories with a sign problem 
and has to be replaced by an alternative mechanism that is worked out in 
detail for QCD in one dimension at nonzero chemical potential. 

PACS numbers: ll.30.Rd, 12.38.Gc 

1. Introduction 

Despite of tremendous efforts [H El El HJ [5] , the phase diagram of QCD 
in the chemical potential temperature plane is only known for fi = and 
at asymptotically large densities. In between these extreme limits, the only 
firm result at low temperatures is the transition to nuclear matter when 
the sum of the chemical potential and the binding energy of nuclear matter 
is equal to the nucleon mass. Most results at intermediate values of the 
chemical potential are based on models such as the Nambu-Jona-Lassinio 
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model [Smj, strong coupling expansions [8l[9], ADS-CFT dualities [13 El 
[T2] and random matrix theory [13^ 111] . 

The reason for the lack of first principle calculations is the phase of 
the fermion determinant at nonzero chemical potential which invalidates 
probabilistic methods. When the sign problem is mild, though, simulations 
can be performed by absorbing the phase in the observable which is known 
as reweighting. The real merit of the work by Fodor and Katz [15] is the 
realization that the sign problem is much less severe close to Tc so that 
sophisticated reweighting methods have a chance to work. In fact, results 
close to Tc and small chemical potential obtained with a variety of methods 
such as imaginary chemical potential [16t I17j . Taylor expansion |18[ I19j . 
reweighting [15] , the canonical ensemble [20] and the density of state method 
|21[ I22j . are in close agreement. 

In this lecture we will discuss the average phase factor of the fermion 
determinant. It provides a direct measure for the severity of the sign prob- 
lem. We will focus on the average phase factor in the microscopic domain 
of QCD where nonperturbative analytical results can be obtained by ex- 
ploiting the equivalence of QCD and chiral random matrix theory [23^ i24j. 
Perturbatively, the average phase factor can be calculated by means of chiral 
perturbation theory up to temperatures close to Tc and chemical potentials 
up to the pion mass [25] . We will find that the sign problem is necessarily 
severe when > m,r/2. The reason is that exponentially large contributions 
to the partition function have to be canceled in order to obtain a chiral con- 
densate that has sensible physical properties such as a discontinuity when 
the quark mass crosses the imaginary axis. 

One disturbing observation is that the relation between the spectral 
density of the Dirac operator and the discontinuity of the chiral conden- 
sate seems to be violated for QCD at nonzero chemical potential. From an 
analysis in the microscopic domain of QCD it was found [26] that the dis- 
continuity arises due to an alternative mechanism. Lending support to its 
universality, the same mechanism is at work [27j for QCD in one dimension. 
This case, because of its simplicity, will be discussed in detail below. 

In this lecture we will show that chiral random matrix theory has added 
significantly to our understanding of chiral symmetry breaking and the sign 
problem for QCD at nonzero chemical potential. This adds to a long list of 
successes of random matrix theory in this field such as the understanding 
of quenched approximation [28], the critical endpoint in QCD [13\ I14j. the 
macroscopic spectral density pS', [3D1 [311 1321 1331 IM], the microscopic 
spectral density [33 [Ml 1311 138 , 39, 401 SI] and Yang-Lee zeros [Ml IMl I13]. 
For more successes we refer to reviews of this subject |44[ I45| . 

After some introductory remarks on QCD at nonzero chemical potential 
in section 2, we will discuss the sign problem in section 3. The microscopic 
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domain of QCD will be introduced in section 4. Results for the average 
phase factor will be presented in section 5, and their relation with the Dirac 
spectrum is examined in section 6. Section 7 contains a detailed discussion 
of the alternative to the Banks-Casher formula for the example of one- 
dimensional QCD at nonzero chemical potential. Concluding remarks are 
made in section 8. 

2. QCD at Nonzero Chemical Potential 

The QCD partition at temperature and chemical potential /x is given 

by 

k 

where the sum is over all states. This partition function can be rewritten 
as a Euclidean quantum field theory 

Zqcd = (ndet(i:' + m/ + ^7o))ym, (2) 
/ 

where the average is over the Yang-Mills action. The Dirac operator is 
denoted by D and the product is over Nj flavors with mass mj. The 
Dirac operator is nonhermitean whereas /i7o is Hermitean so that the total 
Dirac operator has no hermiticity properties. In lattice QCD, the chemical 
potential enters on each time-like link such that forward hopping is enhanced 
by exp(yu) whereas backward hopping is suppressed by exp(— //) [36]. The 
asymmetry between forward and backward propagation also implies that 
the Dirac operator is nonhermitean for n ^ 0: 

det{D + m + /X70) = | det(L' + m + /U7o)|e*^ = Y[{\k + m). (3) 

k 

If the average phase factor vanishes in the thermodynamic limit, Monte 
Carlo simulations are not possible. This problem is known as the sign 
problem. Nevertheless, we emphasize that it is our aim to understand QCD 
at 7^ starting from first principles. 

Let us now discuss the phase diagram in the chemical potential tempera- 
ture plane. A schematic phase diagram is shown in Fig. [TJ Even though the 
11 = axis is rather well-understood, there is still an ongoing controversy 
about the value of the crossover temperature. Whether its value of 190 MeV 
|47j or 170 MeV [48j or somewhere in between, will definitely be resolved 
by future lattice simulations. At asymptotically large values of the chemi- 
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Fig. 1. A schematic phase diagram of QCD in the temperature chemical potential 
plane 



cal potential, perturbative calculations, show that QCD is superconducting 
|49j . We will not further discuss this region. 

At low temperature and /u < (m^r — BE) /3 with BE the binding energy 
per nucleon of nuclear matter, only the vacuum state contributes to the 
QCD partition function, and its free energy is ^u-independent. At /x = 
{niN — BE)/ 3 a transition to nuclear matter takes place. These are the 
only solid results at intermediate densities. All other results at intermediate 
density are model dependent. We agree with McLerran and Pisarki [50j 
that confining forces play an essential role, and that one better relies on 
models where confinement is manifest. One such model is the Skyrme model 
which is believed to be an accurate description of QCD at large Nc- With 
increasing density this model undergoes a transition [5l] to a qualitatively 
different phase. It was shown in [521153] that the dense phase that minimizes 
the energy is a chiral crystal of -B = ^ objects [53] which has restored chiral 
symmetry. This state is strongly bound with a binding energy of at least 
100 MeV per nucleon over nuclear matter. We expect that it takes similar 
temperatures to melt this phase. Most recently such scenario was advocated 
in [50] for the large Nc phase diagram of QCD. If indeed the phase transition 
to the chiral symmetric phase is of first order, the critical line should end 
in a critical endpoint to allow for a crossover transition at small chemical 
potential. 
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3. Triage of the Sign Problem 

A quantitative measure for the severity of the sign problem is given 
by the ratio of the QCD partition function and the phase quenched QCD 
partition function, 

{det^{D + m + ^70)) _ ^-v{Ft,j=2-F^^) 



(|det(D + m + /i7o)|2) 

Because both the numerator and denominator are physical partition func- 
tions they can be expressed in terms of an extensive free energy. At nonzero 
temperature, the difference -FV/=2 — -^pq is always nonzero so that, in the 
thermodynamic limit, the sign problem becomes prohibitively severe. There 
is no reason, though, to doubt that, despite these cancellations, the large 
volume limit is still smooth. What happens is that, with increasing volume, 
it becomes more and more difficult to generate the QCD partition function 
from the phase quenched ensemble. By factoring the determinant of the 
Dirac operator into its absolute and phase factor, exp(i0), the ratio @ can 
be interpreted as the phase quenched expectation value of exp(2i^), 

(det^(D + m + ^7o)) _ , 2ie\ 

\e /pq. \0) 



(|det(D + m + ^7o)|2) 

The phase quenched QCD partition function can be rewritten as the 
expectation value 

(I det(D + m + M7o)P) = (det(L» + m + fijo) det{D + m- ^70)). (6) 

This means that the two flavors have opposite charge with respect to the 
chemical potential. In other words, can be interpreted as an isospin 
chemical potential. Therefore, in the low-temperature limit, the free energy 
of phase quenched QCD remains constant for fj, < m-,r/2. At = 
a phase transition to a Bose-condensed phase of pions takes place so that 
the free energy becomes /i-dependent for ^ > 771^^/2 [55l |56l [57]. In the 
low temperature limit, the free energy of QCD is /^-independent for < 
{rriN—BE) /?,. This implies that the free energy of QCD and phase quenched 
QCD are different for 771,^/2 < /j, < {rriM — BE)/3 resulting in a severe sign 
problem. 

Below we will analyze the average phase factor in the microscopic domain 
of QCD and for one-dimensional QCD. 

4. Microscopic Domain of QCD and Random Matrix Theory 

The QCD Dirac spectrum can be probed by including additional bosonic 
and fermionic determinants in the partition function with quark masses 
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equal to complex parameters that can be varied independent of the QCD 
quark masses. For definiteness lets us assume that we have Nf original 
quarks and s additional quarks with masses given by 

mi,---mNf and zi,---Zs, (7) 

respectively. For fixed QCD quark masses, it is always possible to choose 
these additional quark masses such that the associated Goldstone bosons 
are much lighter than the Goldstone bosons of the original quarks (the 
case of massless quarks has to be treated separately; we assume that all 
quark masses are nonzero). Using the same arguments as given by Gasser 
and Leutwyler for the e-domain of QCD [58], to leading order in the chiral 
expansion of the z^-quarks, the partition function factorizes as [59j 

Z{mi,--- ,mNf,zi,--- ,Zs) = ZgcDimi,--- ,mNf)Z{zi,--- ,Zs) (8) 

if the Compton wave lengths of the additional bosons containing the quark 
masses Zk are much larger than the size of the box. Using the Gell-Mann- 
Oakes-Renner relation, the microscopic domain of QCD is given by ^59j 

« w 

where S is the chiral condensate, F the pion decay constant and V the 
volume of the box. In solid state physics, this scale for z is known as the 
Thouless energy and was discussed in the context of QCD in [601 ED [62] . 
Of course, we can consider QCD with quark masses nik in the microscopic 
domain ([9]) which is also known as the e-domain of QCD. This domain will 
also be called the microscopic domain of QCD. 

In the domain Q, the z^ dependent part of the partition function is 
given by a unitary matrix integral. At /i 7^ invariance arguments lead to 
the following partition function of the following form [55[ |3^ 

= |dC/Sdet^^e^^[str(M[/t+Mtc/)]-iF2ystr[B,c/][B,c/t]^ (10) 

with quark mass matrix given by M = diag(zi, • • • , Zs) and charge matrix 
equal to B = diag(gi, • • • , g^). The superdeterminant and the supertrace 
are denoted by Sdet and Str, respectively. If we have / fermionic quarks 
and s — f bosonic quarks, the integral is over a supergroup with bosonic 
sector equal to the product of U{f) and the positive definite matrices Gl{s — 
f)/U{s — /). For convergence reasons, the quark mass matrix involving the 
bosonic quarks has to be properly adjusted [35] and additional conjugate 
quarks may have to be introduced [63) . 
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Another representation of the partition function (jlOj) is the large- hmit 
of a random matrix model with the symmetries of Z{zi, ■ ■ ■ , Zg). This ran- 
dom matrix model is obtained by replacing the matrix elements of the Dirac 
operator by Gaussian random numbers [6H [65] . 

m iW + fl \ D/T;tAA - „-NTrW'iW 

ivV^ + /X m 



D=[ -'J - ), P(H^)~e-— (11) 



where W is in general an A x (A + v) matrix with u the topological charge. 
As was shown in [66j , the properties of this theory in the microscopic domain 
are not sensitive to the details of the probability distribution. The reason is 
that for A — > oo with mN fixed, the random matrix model has a mass gap, 
so that it becomes a theory of Goldstone bosons dictated by the pattern of 
spontaneous symmetry breaking with partition function given by (jlOp . 

Philosophically, this is important because of the realization that chaotic 
motion dominates the dynamics of quarks at low energy. Practically, this 
is useful because it enables us to use powerful random matrix techniques to 
calculate physical observables. 

5. The Phase of the Fermion Determinant 

Let us now calculate the average phase factor in the microscopic domain 
of QCD. For simplicity we consider the quenched case so that 

/„2i9\ _ I det(£) + m + ^7o) ^ 

^« ^-^detpt + ,„ + ^^o)^ ^^^^ 

Since the average phase factor is the ratio of two partition functions, its 
value is necessarily real and nonnegative. We evaluate the ratio (|12p in 
the microscopic domain of QCD where only the contribution of the zero 
momentum Goldstone modes has to be taken into account. The partition 
function (jl2p has four different Goldstone modes. If we denote the fermionic 
quark by / and the bosonic quark by 6, they are given by 

//, hh, bf, fb. (13) 

The first two are neutral, but because the charge of conjugate quarks is 
opposite to that of regular quarks, the last two have charge ±2. The masses 
of the corresponding Goldstone bosons are thus given by 

mjj = m^, m5ft = m^, m^y- = + 2/i, mj^ = - 2/i. (14) 

When /i < m^/2 only the vacuum state contributes to the partition function 
so that the free energy is a /i-independent constant and the average phase 
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Fig. 2. The microscopic result of the average phase factor (left) and the average 
phase factor for imaginary chemical potential (right). The points with error bars 
in the right figure are lattice results obtained in [S7] . 



factor is given by the product of the square root of the curvatures of the 
Goldstone modes, 

,^2iey ^ (m^ - 2/i)(m^ + 2//) 

For /X > m7r/2 the massless pion Bose condenses, and the free energy be- 
comes //-dependent so that the average phase factor becomes zero in the 
thermodynamic limit. We thus find 

(e2^^) = 0(m.-2|/.|)fl-^V (16) 



In [23] the average phase factor was calculated by means of the complex 
orthogonal polynomial method of [36\ |68] . The result is given by 

{e^'')Nf=o = 1 - 4.fi^Io{m)Ko{m) (17) 

-T^e / dxxe ^i^o I 7^ ) (-^o(a^)"^/l(rn) - xIi(2;)/o("i)) . 

4/i^ Jm V4/i^/ 

The average phase factor can also be calculated using an imaginary chemical 
potential [67J. However, this gives only the first two terms and misses 
the non-analytic term in the above expression. The non-analytic term is 
essential when fi approaches and cancels the analytic term for fi > 

m^/2. 
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Fig. 3. Scatter plot of Dirac eigenvalues and /i 7^ (left), and in the right panel, 
we show the mass dependence the chiral condensate for ftiU QCD (red curve) and 
quenched or phase quenched QCD (green curve). 



In Figure [2] we show the analytical result for the average phase factor 
weighted with the determinant for two flavors (left figure). The solid curve 
is the mean field result. In the right figure we compare the quenched average 
"phase factor" for imaginary chemical potential with lattice simulations on 
an 8'^ lattice (see [67J). The solid curve is given by the analytical continua- 
tion (/x^ — > — /U^) of the first two terms in ()17p . 

6. Average Phase Factor and Dirac Spectrum 

The sign problem becomes inevitable when the quark mass is inside the 
domain of eigenvalues. This should be obvious because for eigenvalues close 
to the mass, small variations of the gauge field result in large phase fluc- 
tuations. In figure [3] we illustrate the distribution of the Dirac eigenvalues. 
The width of the spectrum can be obtained from chiral perturbation theory 
|34j or chiral random matrix theory [28j . In the quenched case or the phase 
quenched case the chiral condensate at m can be interpreted as the planar 
electric field at m of charges located at the position of the eigenvalues. El- 
ementary electrostatics dictates that the chiral condensate behaves as the 
green curve in the right panel of Figure [3l In particular, there is no discon- 
tinuity at m = 0. Because the low temperature limit of the free energy of 
full QCD does not depend on the chemical potential, the chiral condensate 
of full QCD should have a discontinuity at m = for /x < m^v/S. What has 
to happen is that the phase of the fermion determinant has to cancel the 
decrease in free energy that takes place in the quenched or phase quenched 
theory for [x > m,r/2. Therefore, we necessarily have a severe sign problem 
in this domain. 

The puzzle that the chiral condensate remains the same as the result 
of strong cancellations is known as the "Silver Blaze Problem" [69j. This 
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name has been inspired by the title of a Sherlock Holmes novel by Conrad 
Doyle. The analogy is that the dog did not bark although the racing horse 
"Silver Blaze" disappeared. 

In the literature it has been stated that average spectral density of full 
QCD will show an accumulation of spectral density on the imaginary axis 
consistent with the Banks-Casher formula. We now understand that this is 
not the case. The average spectral density for full QCD was evaluated in the 
microscopic domain of full QCD using the method of complex orthogonal 
polynomials [36]. It was found that it has oscillations in a macroscopic 
region of the complex plane with an amplitude that increases exponentially 
with the volume and a period that is inversely proportional to the volume 
|38j . The cancellations of the exponential large contributions result in a 
condensate that has a discontinuity in the thermodynamic limit [26j. 

The same phenomenon occurs for QCD in one dimension. Although the 
eigenvalues of the Dirac operator are located on an ellipse instead of being 
scattered in the complex plane, the mechanism for generating a discontinuity 
in the chiral condensate is the same as for QCD. To simplify the argument 
we will restrict ourselves to giving a detailed discussion for the case of QCD 
in one dimension only. 



7. One-dimensional QCD 



All gauge fields in one-dimensional QCD can be gauged away with the 
exception of the gauge field at the boundary of the manifold. Therefore 
one-dimensional lattice QCD is given by the random matrix theory 



dUdet^fD 



(18) 



with Dirac operator for n lattice points given by. 



/ ml 



D 



ml 



V 





-efU/2 



ml 







ml / 



(19) 



The chemical potential can also be gauged to the boundary. The full theory 
does not have charged excitations resulting in a chiral condensate that is 
//-independent. The phase quenched theory, on the other hand, undergoes 
a phase transition at /u = m-,j/2 resulting in a different free energy. This 
can be easily shown in the limit of large n by using the explicit result for 
the determinant of the Dirac operator (jl9p 



det D = 2~ 



' det e 



+ e 



-n^c ^gn/^c ^gnMc^^g- 



U^. (20) 
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Fig. 4. The eigenvalue distribution of the one-dimensional QCD Dirac operator 
at nonzero chemical potential and large (left). The graphs show the mass 
dependence of chiral condensate of the partially quenched theory (middle) and the 
full theory (right). 



Here, we introduced the critical chemical potential //c given by the relation 
sinh/ic = fn. For gauge group U{N) the integrals are particularly simple, 
and we obtain in the limit of large n 

Fni=2 = -nNf\^c\, FNf=2.pq = -nNf\fic\ - nNf{\fi\ - \fic\)0{\fi\ - |/ic|), 

(21) 

resulting in a chiral condensate with a mass dependence as shown in Figure 
m Although nothing happens to the Nf = 2 free energy, this comes only as 
a result of exponentially large cancellations in the partition function. This 
is the "Silver Blaze Problem" [BDj mentioned before. We will now illustrate 
how this problem manifests itself in spectrum of the Dirac operator and the 
chiral condensate. 



For simplicity we will only consider the case A'^^ = 1 with U{1) as gauge 
group. Because the phase angles of the eigenvalues are uniformly distributed 
along an ellipse with semi-minor axis equal to sinh fi and semi-major axis 
equal to cosh// (See Figure HI Notice that the figure is for large Nc- For 
{7(1) the phase angles of the eigenvalues are equally spaced.), the eigenvalue 
density of the Nf = 1 theory is given by 
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with the original variables z parameterized as 

z = i(e'^+*" - e-^'-^°), r > 0, a G [0, 27r]. (23) 
The chiral condensate for Nf = 1 given by 



(rifclAfc + m)) 



can then be expressed as 



^Nf=l — / « z ; 

z + m 

da e^^" + e""'^'' — e"^*""'"'^) — g-"(*"+M) 
27r (e"'/^'= + e-"/^'=)(m + (e/^+*° - e-/^-*°^)/2) 

tanh(n^c) 
cosh /ic 



(25) 



When |/i| < |//c|, all four terms in the numerator contribute to the integrand. 
When > \nc\, i-e. the domain where the quark mass is inside the ellipse 
of eigenvalues, the quenched chiral condensate, given by the first two terms 
of the numerator, is zero resulting in a chiral condensate that behaves as 
the phase quenched theory (see middle panel of Figure [H Therefore the 
condensate for |//| > \fic\ is due to the oscillating terms in the numerator. 
A finite result is obtained because the oscillations cancel the exponential 
growth with n of the amplitude. In the thermodynamic limit the tanh /ic 
results in a discontinuity at m = (see Figure H]) . 



8. Conclusions 

Random matrix theory has been invaluable for understanding QCD at 
nonzero chemical potential. In this lecture we have discussed applications 
involving the phase of the fermion determinant at nonzero chemical poten- 
tial. We have shown that QCD has a severe sign problem if the quark mass 
is inside the domain of eigenvalues. In this domain strong cancellations 
lead to a free energy that does not depend on the chemical potential and a 
chiral condensate that has a discontinuity when the quark mass crosses the 
imaginary axis. The latter happens without an accumulation of eigenvalues 
on the imaginary axis, but due to oscillations in the spectral density with 
an amplitude that increases exponentially with the volume. This mecha- 
nism occurs both in the random matrix limit of QCD and in one dimen- 
sional QCD which strongly suggests that it is the generic replacement of 
the Banks-Casher formula for theories with a sign problem. 
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